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ABSTRACT 

The  construction  of  very  good  hyperspectral  sensors  operating  in  the  thermal  infrared  bands  from  8  to 
12  microns  arouses  much  interest  for  the  development  of  data  exploitation  tools.  Temperature  emissivity 
separation  (TES)  algorithms  are  very  important  components  of  a  future  toolbox,  because  they  make  it  possible 
to  extract  these  two  fundamental  targets'  parameters.  The  emissivity  relies  on  the  nature  of  the  target's  surface 
materials,  while  the  temperature  gives  information  related  to  their  use  and  relationship  with  the  environment. 
The  TES  technique  presented  in  this  paper  is  based  on  iteration  on  temperature  principle,  where  a  total 
square  error  criterion  is  used  to  estimate  the  temperature.  The  complete  procedure  is  described  in  the  paper. 
Its  sensitivity  to  noise  is  studied  and  a  mathematical  behavior  model  is  provided.  The  model  is  validated 
through  a  Monte-Carlo  simulation  of  the  technique's  operation. 


1.0  INTRODUCTION 

In  the  long  wavelength  infrared  (LWIR)  band,  extending  from  8  to  12  microns,  the  temperature  and  the 
emissivity  can  be  defined  as  the  fundamental  parameters  of  the  imaged  material  or  targets.  The  emissivity 
provides  information  related  to  the  target's  nature  while  the  temperature  relates  to  its  relationship  with  the 
environment  or  to  its  activity.  Another  very  interesting  feature  of  the  LWIR  bands  is  that  the  imagery  could  be 
acquired  by  day  or  night. 

The  processing  chain  surrounding  calibrated  thermal  infrared  hyperspectral  imagery  leading  to  the  extraction 
of  the  imaged  material  fundamental  properties  are  atmospheric  compensation  followed  by  temperature- 
emissivity  separation.  The  atmospheric  compensation  is  the  process  by  which  the  atmospheric  transmittance 
and  path  radiance  are  removed  from  the  imagery.  This  step  provides  two  important  results:  the  ground-leaving 
radiance  and  the  atmospheric  downwelling  irradiance.  The  path  radiance  is  the  energy  generated  by  the 
atmosphere  on  the  path  from  the  target  to  the  sensor.  The  transmittance  is  the  amount  of  energy  emitted  by  the 
ground  and  lost  on  the  path  from  the  target  to  the  sensor.  The  downwelling  irradiance  is  the  energy  incident  on 
the  target  originating  from  the  hemisphere  above  the  target;  it  is  often  converted  to  radiance  by  assuming 
lambertian  reflection  from  the  target.  Finally,  the  ground-leaving  radiance  is  the  energy  leaving  the  target  and 
is  the  combination  of  the  target's  self  radiation  and  reflection  of  the  atmospheric  downwelling  irradiance. 


Lahaie,  P.  (2005)  A  Temperature  and  Emissivity  Separation  Technique  for  Thermal  Hyperspectral  Imagers.  In  Emerging  EO 
Phenomenology  (pp.  4-1  -  4-14).  Meeting  Proceedings  RTO-MP-SET-094,  Paper  4.  Neuilly-sur-Seine,  France:  RTO.  Available 
from:  http://www.rto.nato.int/abstracts.asp. 


RTO-MP-SET-094 


4-1 


UNCLASSIFIED/UNLIMITED 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

01  OCT  2005 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

A  Temperature  and  Emissivity  Separation  Technique  for  Thermal 
Hyperspectral  Imagers 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROIECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Data  Exploitation  Group  Optronic  Surveillance  Section,  DRDC 

Valcartier  2459  Pie-XI  Blvd.  North  Val-Belair,  Quebec 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM202114,  RTO-MP-SET-094.  Emerging  EO  Phenomenology  (Naissance  de  la  phenomenologie 
et  de  la  technologie  electro-optique).,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBIECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

30 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


UNCLASSIFIED/UNLIMITED 


A  Temperature  and  Emissivity  Separation 
Technique  for  Thermal  Hyperspectral  Imagers 


ORGANIZATION 


Using  the  downwelling  irradiance  and  the  ground-leaving  radiance,  it  is  possible  to  extract  the  emissivity  and 
the  temperature  of  a  target.  The  downwelling  irradiance  acts  in  the  process  as  the  reference.  This  principle  has 
been  used  by  Borel  [1]  and  resulted  in  the  ISSTES  algorithm.  This  algorithm  has  subsequently  been  studied 
thoroughly  by  Ingram  and  Muse  [2],  In  our  technique,  we  use  the  downwelling  irradiance  as  the  reference  for 
the  emissivity.  We  also  use  iteration  for  the  determination  of  temperature,  like  what  is  done  in  Borel's 
procedure.  Our  technique's  difference  from  ISSTES  lies  in  the  method  used  for  selecting  the  right  temperature 
and  its  corresponding  emissivity.  That  difference  leads  to  an  increase  in  resistance  to  noise  and  to  impairments 
such  as  the  wrong  estimation  of  the  downwelling  irradiance. 

The  paper  is  arranged  as  follows.  In  section  2  the  assumed  signal  is  described.  Section  3  describes  the  general 
algorithm  and  highlights  its  most  important  features.  Section  4  is  devoted  to  a  sensitivity  analysis  for  noise 
and  section  5  provides  some  experimental  results  obtained  from  ground  sensing.  Finally,  the  conclusion  is 
presented  in  section  6. 

2.0  SIGNAL  MODEL 

A  complete  signal  model,  taking  into  account  every  features  of  the  atmosphere,  is  very  complicated.  It 
involves  phenomenon  such  as  heat-generated  signal  by  ground  objects  and  the  atmosphere,  scattering  by  the 
atmosphere  (aerosols  and  clouds),  absorption  by  the  surface  and  by  the  atmosphere.  If  every  single 
phenomenon  was  taken  into  account  perfectly,  a  forward  computation  could  be  possible,  but  the  inversion 
rapidly  end  into  an  intractable  problem.  One  solution  to  that  problem  is  to  specify  conditions  for  which  the 
model  possesses  a  simpler  solution.  A  more  tractable,  but  still  very  complicated  situation  arises  when  a  clear 
sky  is  the  weather  condition.  For  a  given  spectral  band,  the  signal  model  becomes: 
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where  a  is  the  wavenumber  scale.  This  scale  will  be  used  throughout  this  paper  since  we  also  use  it  in 
conjunction  with  MODTRAN  and  with  the  instruments  used  to  acquire  data  (ABB-Bomem  Michelson 
interferometer).  The  term  Lmn  is  the  measured  radiance  in  band  n  and  fn  is  a  normalized  weighting  function  for 
the  spectral  response  of  the  instrument,  s  is  the  emissivity  of  the  target,  BT  is  the  Planck  function  at 
temperature  T.  Li  is  the  atmospheric  radiance  incident  on  the  surface  from  the  sky  from  direction  {0,  (j),  r  is 
the  atmospheric  transmittance  on  the  path  from  the  target  to  the  sensor  and  Lp  is  the  path  radiance  accumulated 
from  the  target  to  the  sensor  and  is  given  by: 


Lp=\B(T(y),a 


exp 
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where  T(y)  is  the  temperature  of  the  atmosphere  at  elevation  v  and  a  is  the  extinction  coefficient  of  the 
atmosphere  at  elevation  y.  In  this  model  the  influence  of  single  scattering  by  aerosol  particles  is  neglected  as  a 
contribution  to  the  signal.  The  aerosols'  effects  are  considered  however,  in  the  extinction  coefficient.  Another 
simplification  takes  place  in  the  reflective  characteristics  of  the  targets  that  are  considered  to  be  lambertian. 
Directional  reflection  peculiarities  of  targets  are  discarded  in  the  model.  The  mathematical  simulation  of  a 
potential  sensor  measurement  therefore  requires:  the  knowledge  of  the  sensor's  characteristics  such  as  the 
spectral  response  of  each  of  its  bands;  the  target's  characteristics,  which  are  the  spectral  emissivity,  reflection 
characteristics  and  its  temperature;  and  the  complete  knowledge  of  the  atmosphere,  which  are  the  temperature, 
the  pressure  and  the  humidity  gradients. 
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The  clear  sky  model  needs  to  be  simplified  even  more  with  the  use  of  some  physical  assumption.  The  targets' 
emissivities  can  be  assumed  constant  within  the  bounds  of  each  band.  This  holds  also  for  the  blackbody 
function.  This  means  that  a  bulk  transmittance  can  be  used  for  the  whole  band.  The  downwelling  irradiance  is 
spectrally  highly  variable.  The  bulk  transmittance  for  the  computation  of  its  reflection  transmission  through 
the  atmosphere  cannot  be  used  if  high  precision  is  required.  However,  if  the  sky  is  clear  and  the  atmosphere 
equivalent  temperature  is  much  smaller  than  the  ground  temperature,  in  conjunction  with  generally  high  target 
emissivity  renders  marginal  the  contribution  of  downwelling  irradiance.  The  bulk  transmittance  can  therefore 
be  used  as  an  approximation  for  the  reflected  downwelling  irradiance  contribution.  This  will  introduce  more 
difficulties  for  processing  of  highly  reflective  targets,  but  since  these  targets  are  already  extremely  difficult  to 
process,  it  does  not  increase  the  burden  too  heavily.  The  path  radiance  is  additive  and  can  therefore  be 
considered  separately.  For  any  given  band,  the  model  can  be  simplified  to  provide  the  following  equation: 

r  t  \ 


e.B.(T)  +  ( 1-0- 


+  Lpn+N 
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In  the  remainder  of  the  document,  the  n  factor  dividing  the  Ls  contribution  of  atmospheric  downwelling 
irradiance  will  be  omitted.  It  must  however  be  kept  in  mind  that  the  Ls  variable  is  been  converted  into  radiance 
through  this  transformation.  The  bracketed  part  of  the  preceding  equation  is  known  as  the  ground-leaving 
radiance.  The  N  variable  is  the  sensor  added  noise  for  the  band  of  interest.  The  removal  of  the  path  radiance 
and  of  the  transmittance  constitutes  the  first  component  of  the  processing  chain  i.e.  the  atmospheric 
compensation.  This  provides  the  ground-leaving  radiance.  This  last  quantity  can  be  measured  directly  by  a 
ground  spectrometer  and  this  way  data  can  be  gathered  to  verify  specifically  a  temperature-emissivity 
separation  algorithm.  The  signal  model  used  in  the  remainder  of  the  paper  is  therefore: 

L„=e,B,(T)  +  (  1 (4) 

3.0  ALGORITHM  DESCRIPTION 

3.1  Introduction 

The  TES  system  is  built  around  the  minimization  of  an  error  function.  The  error  function  can  be  the  total 
square  error  or  the  total  absolute  error.  In  the  total  square  error  case,  an  error  is  computed  for  each  band  and 
the  squares  of  band  results  are  added  together.  For  the  total  absolute  error,  the  absolute  values  of  the  error  in 
each  band  are  added  together.  The  system  is  assumed  linear  in  the  region  of  the  minimum  error  and  for  this 
reason  the  two  error  computation  methods  should  be  unbiased  and  therefore  generate  similar  minimum  values. 
However,  there  could  be  a  difference  in  performance  especially  for  the  variance  of  the  error  on  temperature 
due  to  noise.  We  prefer  the  total  square  error  because  in  that  case  it  is  possible  to  get  an  expression  for  the 
computation  of  the  temperature  error  variance.  An  analytical  expression  is  not  possible  for  the  total  absolute 
value,  and  any  verification  must  be  done  using  a  Monte-Carlo  simulation.  This  difficulty  makes  the 
elaboration  of  a  design  procedure  for  components  of  the  technique  very  difficult.  The  total  square  error  is 
given  by: 

n 

where  Rgn  is  the  ground-emitted  radiance  in  the  «th  band  and  Rf„  is  the  emissivity-filtered  radiance  computed 
for  band  n.  The  algorithm  structure  for  a  given  pixel  is  given  in  figure  1.  It  assumes  that  the  ground-leaving 
radiance  and  the  downwelling  irradiance  have  been  provided  as  inputs.  Prior  to  the  operation  on  image's  pixels 
the  algorithm  initialization  must  be  done.  This  step  encompasses  the  computation  of  pixels'  ground-leaving 
radiance  and  the  downwelling  irradiance  for  the  image.  Figure  1  shows  the  structure  of  the  algorithm  for 
pixel-by-pixel  processing.  It  contains  the  following  steps:  computation  of  a  start  temperature;  emissivity 
computation;  emissivity  smoothing;  computation  of  radiance,  using  smoothed  emissivity;  computation  of  the 
total  squared  error  on  radiance;  decision  to  stop  the  process  and  finally  estimation  of  a  new  trial  temperature. 
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Figure  1:  Structure  of  the  proposed  algorithm 


3.2  Algorithm  initialization 

This  step  is  not  included  in  figure  1 .  The  initialisation  phase  of  the  algorithm  is  constituted  by  the  input  of 
basic  data,  which  are  the  image  size,  the  atmospheric  profiles  or  atmospheric  optical  parameters 
(transmittance,  path  radiance  and  downwelling  irradiance),  the  sensor  characteristics  such  as  the  number  of 
bands,  their  centers  and  width  and  possibly  there  shift  with  lateral  position  also  known  as  the  smile  and  the 
sensor  altitude.  If  the  atmospheric  profiles  were  provided,  it  includes  the  computated  atmospheric  optical 
parameters.  The  computation  is  skipped  if  the  optical  atmospheric  parameters  for  each  band  are  provided 
instead  of  the  atmospheric  profiles.  The  image  atmospheric  compensation  that  comprises  the  elimination  from 
the  data  of  the  path  radiance  and  of  the  transmittance  is  part  of  the  algorithm  initialisation.  Another  future 
component  of  the  initialization  is  the  computation  of  the  filter's  coefficients  according  to  the  weather  and 
geographical  location.  This  step  could  be  added  in  future  when  a  filter  design  technique  will  exist. 
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3.3  Pixel  initialisation 

Figure  1  illustrates  the  algorithm  work  for  pixel  processing  that  is  the  object  of  TES.  The  pixel  processing 
needs  to  be  started  at  the  best  possible  temperature  estimation.  The  constraint  associated  with  the  rough 
estimate  of  the  temperature  is  that  it  must  be  within  the  region  of  stability  of  the  algorithm  and  thus  the 
estimate  must  be  as  close  as  possible  to  the  true  temperature  of  the  pixel's  sample.  The  technique  uses  two 
adjacent  or  nearly  adjacent  bands  of  ground-leaving  radiance  and  their  corresponding  downwelling  irradiance. 
The  assumptions  are:  the  noise  can  be  neglected  and  the  emissivities  of  the  sample  in  the  two  bands  are  so 
similar  that  they  can  be  assumed  equal.  The  proximity  of  the  two  bands  justifies  this  approximation.  One  very 
important  characteristic  for  the  selection  of  the  two  bands  is  that  there  should  be  a  high  difference  in  radiance 
due  mainly  to  the  difference  in  downwelling  irradiance.  The  assumption  that  the  emissivities  are  the  same 
reduces  the  number  of  unknowns  to  two  and  therefore  the  problem  could  be  solved  with  some  small  error.  The 
two  equations  become: 


=sBx{T)  +  {\-s)Lx 

R , 


(?) 


•2  =  sB2  (T)  +  (1  -  s)L2 

where  the  index  of  the  blackbody  functions  indicates  one  band  and  it's  nearest  neighbour.  Another 
approximation  can  be  done  that  assumes  the  blackbody  functions  for  the  two  bands  are  equal.  The  equation  set 
becomes  linear  and  the  estimate  for  temperature  becomes: 

f  nr  nr  \ 
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3.4  Emissivity  computation 

The  emissivity  computation  is  the  first  step  performed  in  the  iteration  loop  on  temperature,  as  given  by 
equation  (9). 


Rs-L 

B-L 


(9) 


3.5  Emissivity  filtering 

Emissivity  filtering  is  one  of  the  critical  parts  of  the  algorithm.  The  filter  behaviour  determines  to  a  large 
extent  the  capability  of  the  algorithm  to  estimate  the  temperature  with  precision.  In  previous  versions  of  the 
algorithm  the  filter  was  constituted  around  a  smoothing  process.  The  fitting  of  a  polynomial  and  the 
computation  results  of  a  smoothed  version  of  the  emissivity  were  used  afterwards  to  compute  the  error.  This 
technique  is  reasonably  good  for  the  case  of  emissivities  containing  only  small  variations.  If  an  emissivity 
happens  to  be  very  variable  the  order  of  the  polynomial  required  to  fit  it  becomes  too  large  so  the  coefficients 
cannot  be  estimated  accurately.  The  method  is  bounded  by  the  use  of  a  polynomial  order  lower  or  equal  to  5. 

In  the  analysis  of  that  method  it  has  been  observed  that  the  smoothing  process  is  linear  and  can  be  described 
by  the  use  of  a  linear  filter  applied  to  the  emissivity  vectors.  This  suggested  the  use  of  such  filter  for 
application  on  the  emissivity.  However,  neither  a  theory  for  the  result's  evaluation  or  the  design  of  such  filters 
has  been  identified  in  the  literature  in  this  kind  of  application.  The  application  of  the  filter  on  the  emissivity  is 
given  by: 

s  =  Gs  (10) 

where  s  is  the  emissivity,  G  is  the  matrix  filter,  and  £  is  the  filtered  emissivity. 


3.6  Radiance  estimation 

The  radiance  estimation  is  a  simple  step  performed  once  the  emissivity  has  been  filtered.  The  equation  is  the 
following: 
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Rf  =  sB  +  (l  -  s]l  (1 1) 

where  Rf  is  the  radiance  computed  with  the  filtered  emissivity,  B  is  the  blackbody  function  computed  at  the 
trial  temperature  and  £  is  the  filtered  emissivity. 

3.7  Error  estimation 

The  error  estimation  is  calculated  using  the  following  expression  for  the  total  square  error. 

e2=ZK-rJ  <i2> 

It 

If  the  total  absolute  error  is  preferred  the  equation  is: 

£  =  2X-*/.|  <13> 

n 

3.8  Decision 

The  decision  to  terminate  the  pixel  processing  or  to  modify  the  temperature  is  based  on  the  simplex  method. 
The  algorithm  is  initialized  with  a  step  on  temperature  that  is  set  at  one  degree.  Initially  the  temperature  is 
supposed  to  be  a  small  amount  under  the  computed  minimal  temperature  so  the  initial  direction  of  movement 
for  the  temperature  is  upward.  Each  time  a  temperature  is  updated,  the  error  must  be  computed.  If  the  error 
decreases,  the  temperature  direction  and  step  are  maintained.  When  the  error  increases,  even  just  a  little  bit, 
the  temperature  movement  direction  is  reversed  and  the  step  is  decreased.  That  procedure  is  followed  until  the 
temperature  step  reaches  a  minimum  value.  The  temperature  generating  the  minimum  error  value  at  the 
smaller  temperature  step  is  selected  as  the  pixel  temperature.  Figure  2  shows  graphically  this  procedure. 


Figure  2:  Description  of  the  minimum  search  algorithm. 
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The  algorithm  operation  has  been  described  as  a  system.  Improvements  could  probably  be  added  in  the 
general  operation  such  as  the  variation  in  temperature  steps  or  the  temperature  adaptation  could  be  modified 
using  a  different  procedure  with  the  use  of  Newton-Raphson  method  or  steepest  descent  algorithm.  In  any 
case  this  constitutes  a  rather  minor  modification  compared  to  the  design  of  adequate  filters  or  the  use  of  a 
different  criterion  for  temperature  selection. 

4.0  NOISE  SENSITIVITY  ANALYSIS 

4.1  Introduction 

Sensor  noise  is  the  only  impairment  that  will  always  be  present  in  any  measurement  system.  All  other 
problems  could  be  avoided  or  greatly  reduced  by  modeling  and  processing.  These  other  problems  are  related 
to  the  estimation  of  the  required  variables  of  the  system  (the  downwelling  irradiance,  the  transmittance,  the 
path  radiance  and  the  emissivity  variations).  It  is  reasonable  to  think  that  the  optical  atmospheric  parameters 
as  well  as  the  sensor's  parameters  will  be  thoroughly  known  for  a  given  sensor  and  measurement  setup  so, 
here,  we  limit  ourselves  to  the  evaluation  of  the  impacts  of  sensor  noise  on  the  temperature  and  emissivity 
estimation.  What  is  required  are  the  temperature  introduced  bias  and  the  variance  of  the  temperature 
estimation  due  to  noise. 


4.2  Temperature  bias  due  to  noise 

We  begin  with  the  total  square  error  as  a  function  of  temperature. 
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Assuming  the  blackbody  function  can  be  expanded  in  a  Taylor  series  for  the  region  near  the  temperature  of 
interest,  the  total  square  error  containing  noise  can  be  expressed  by: 
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Using  the  fact  that  B(T )  «  Bo  +  B]Td 

Where  B0  is  the  blackbody  function  at  the  temperature  of  the  target  and  5/  is  the  first  order  derivative  of  the 
blackbody  function  estimated  at  the  target's  temperature.  Td  is  the  deviation  from  that  temperature  that  is 
introduced  by  noise  and  sQ  is  the  true  emissivity  of  the  target. 


This  function  is  quadratic  on  temperature  and  is  valid  in  the  region  around  the  target's  temperature.  In  this 
function  all  the  variables  are  fixed  or  supposed  to  be  known  in  a  given  case.  Then  to  find  the  best  estimate  for 
the  temperature  one  only  has  to  minimize  the  function  by  differentiating  the  error  function  with  respect  to 
temperature  and  equating  the  derivative  to  zero  to  extract  the  temperature  as  a  function  of  all  the  other 
variables.  With  a  little  rearrangement  the  derivative  becomes. 
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This  expression  is  very  complicated  and  a  lot  of  terms  contained  inside  it  possess  a  low  magnitude  compared 
to  major  terms.  A  study  of  each  term's  behaviour  and  magnitude  near  the  minimum  temperature  is  needed. 
The  contrast  between  the  ground  and  the  sky  radiances  is  the  signal  component  possessing  the  highest 
magnitude  in  the  above  expression.  Following  that  magnitude,  are  second  order  components  such  as  the  noise 
and  the  radiance  difference  due  to  the  difference  in  temperature  Td.  The  emissivity  possesses  a  magnitude  of  1 
and  the  temperature  difference  is  expected  to  be  very  small  if  B]Td  is  considered  to  have  a  magnitude  at  least 
comparable  to  the  noise.  Neglecting  the  smaller  terms  with  the  use  of  the  preceding  arguments  and  inverting 
equation  17  for  temperature  difference  leads  to: 
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This  expression  provides  the  bias  due  to  noise  in  a  given  single  case,  where  the  noise  would  be  known.  The 
temperature  estimate  is  unbiased  if  the  noise  magnitude  is  low  enough  and  if  the  noise  mean  is  null.  To  be 
valid  expression  18  requires  a  high  signal  to  ratio. 


4.3  Temperature  Variance  due  to  noise 

The  variance  of  the  temperature  estimation  can  be  estimated  directly  from  these  equations  and  gives: 
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4.4  Results  comparisons  of  simulations  and  computations 


Noise  equivalent  radiance  [W/(m2.sr.cm'1)] 


Figure  3:  Example  of  a  noisy  signal  in  the  first  graph  for  which  the  red  curve  radiance  is  computed 
using  a  unitary  emissivity  and  noise  of  a  high  magnitude  is  added.  Second  graph  shows  the  bias  on 
the  estimated  temperature  for  a  case  where  the  noise  standard  deviation  is  constant  and  indicated  on 
the  abscissa,  is  the  same  for  each  band  and  the  basis  signal  is  computed  using  a  0.96  constant 
emissivity  and  a  temperature  of  300K.  The  downweiiing  irradiance  is  computed  using  MODTRAN 

standard  tropical  atmosphere. 


Figure  4:  The  left  graphics  shows  a  comparison  between  the  standard  deviation  computed  with 
equation  19  and  a  MODTRAN  standard  tropical  atmosphere  and  a  Mont-Carlo  simulation  in  the  same 
condition.  The  temperature  of  the  sample  is  300K.  Right  graphics  shows  the  impact  of  emissivity 

value  on  the  temperature  standard  deviation. 


Figures  3  and  4  show  some  of  the  performance  of  the  algorithm  in  simulations.  In  these  simulations  and 
estimation,  there  is  no  consideration  of  the  effect  on  the  estimated  temperature  for  errors  on  atmospheric 
optical  parameters.  The  temperature  of  the  samples  is  300K  and  the  atmosphere  is  the  MODTRAN  standard 
tropical  atmosphere.  The  contrast  is  therefore  very  low  for  these  conditions  since  it  is  expected  that  for 
daylight  condition  the  temperatures  of  surface  materials  should  be  much  higher.  The  tropical  atmosphere  is 
also  very  hot  and  wet.  It  is  believed  according  to  Ingram  and  Muse  that  this  represents  challenging  conditions 
for  the  algorithm.  This  is  in  agreement  with  the  downweiiing  irradiance  of  the  left  hand  graphics  of  figure  3 
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where  the  contrast  is  effectively  very  small.  The  left  graph  of  figure  4  and  the  right  graph  of  figure  3  shows 
that  expression  19  could  effectively  be  used  and  yield  good  results  for  the  estimation  of  the  temperature 
estimate  standard  deviation  once  the  atmospheric  conditions  and  sensor's  characteristics  are  well  known  and 
when  the  signal  to  noise  ratio  is  good. 


Figure  5:  Left  graph  represents  the  impact  of  the  filter's  characteristics  (filter  width)  on  the  estimation 
of  the  temperature.  The  emissivity  is  constant  at  0.96.  The  right  graph  is  the  standard  deviation  for  the 
emissivity  given  a  mean  signal  to  noise  ratio  for  the  whole  sensor  signal. 

The  width  of  a  filter  represents  its  capability  to  follow  fast  variations  in  the  emissivity.  The  narrower  it  is  the 
faster  it  can  adapt  to  fast  variations.  The  estimation  of  the  temperature  for  such  a  fast  varying  emissivity  will 
therefore  be  better  using  a  narrow  filter.  However,  the  right  graph  of  figure  4  also  shows  that  the  wider  a  filter 
is,  the  smaller  the  standard  deviation  will  be.  A  compromise  should  therefore  be  made  relatively  to  the  filter 
that  will  be  used  in  a  given  situation.  It  is  interesting  to  note  that  in  the  formalism,  there  is  no  requirement  for 
all  the  lines  of  a  filter  to  be  equivalent.  It  means  that  some  part  of  the  filter  can  be  faster  to  adapt  while  other 
parts  are  slower.  This  complicates  enormously  the  filter  design  by  providing  a  very  high  level  of  freedom  for  a 
design  procedure.  One  should  therefore  be  very  clever  in  specifying  constraints. 

5.0  GROUND  MEASUREMENT  EXAMPLE 

In  this  section  we  show  some  results  obtained  through  ground  measurements  during  June  of  2005.  These 
measurements  are  mostly  used  to  show  that  the  signal  model,  assuming  lambertian  reflections,  can  be  used  in 
a  general  way.  They  were  performed  in  an  open  area  of  DRDC-Valcartier,  where  there  are  no  high  building  in 
the  direct  vicinity. 

The  measurement  setup  consists  of  an  ABB-Bomem  MR300  spectrometer  mounted  on  a  pan  and  tilt  device. 
The  spectrometer  aims  a  folding  mirror  that  reflects  the  radiance  from  a  target  placed  on  the  ground.  The 
measurements  were  done  on  a  sunny  day  during  the  morning  and  the  afternoon  in  two  different  runs.  To 
calibrate  the  spectrometer  three  measurements  are  made,  to  obtain  the  gain  and  the  offset  of  the  instrument 
and  of  the  various  components.  A  hot  blackbody  source  is  first  measured,  then  the  downwelling  irradiance  is 
measured  with  an  infragold  plate  and  an  ambient  temperature  blackbody  source  is  used.  The  downwelling 
irradiance  measurement  is  used  to  compensate  for  the  small  reflections  generated  by  the  blackbodies,  which 
possess  an  emissivity  near  98.5%. 

Three  measurements  examples  done  during  the  morning  and  the  afternoon  are  shown  at  figure  6  to  8.  The  first 
target  is  a  sandblasted  aluminium  panel  used  in  airborne  trials  for  calibrating  sensors.  The  second  target  is  a 
plain  plywood  panel  representing  typical  construction  material  and  the  third  is  a  roofing  material  panel. 
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Calibrated  radiance 


Calibrated  downwelling  irradiance 


Wavenumber  [cm'1] 
Aluminum  emissivity 


Figure  6:  Aluminium  panel:  Panel  photograph,  ground-leaving  radiance,  downwelling  irradiance  and 

emissivity.  Estimated  temperature:  AM  65°C,  PM  59°C 


Calibrated  radiance 


Calibrated  downwelling  irradiance 


Wavenumber  [cm'1] 
Plywood  emissivity 
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Figure  7:  Plywood  Panel:  Photograph,  ground-leaving  radiance,  downwelling  irradiance  and 
emissivity.  Estimated  temperature:  AM  38.1°C  PM  35.7°C 
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Calibrated  radiance 
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Calibrated  downwelling  irradiance 
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Figure  8:  Roofing  material:  Photograph,  ground-leaving  radiance,  downwelling  irradiance  and 
emissivity.  Estimated  temperature:  AM  61.5°C,  PM  54.5°C 

The  temperature  estimation  given  in  the  figure  caption  is  the  temperature  at  which  the  emissivity  has  been 
computed  for  each  samples.  Temperature  measurements  for  validation  are  notoriously  difficult  to  obtain, 
using  any  kind  of  technique.  For  instance  a  radiometer  would  not  provide  accurate  measurement  for  a  sample 
such  as  the  roofing  material  since  the  emissivity  is  not  constant.  A  thermocouple  placed  at  the  surface  may 
find  a  different  temperature  to  other  parts  of  the  surface  since  conductivity  and  cooling  by  convection  due  to 
wind  may  be  different.  A  measurement  done  underneath  the  material  could  also  show  a  large  difference  with 
the  surface  because  of  the  time  lag  due  to  the  thickness  of  the  material.  The  contact  between  the  material  and 
the  substrate  could  also  interfere  with  the  measurement.  These  various  problems  complicate  the  task  of 
validating  the  algorithms  using  independent  measurements.  The  correspondence  of  the  emissivity  for  each 
measurement  is  clearly  seen  from  the  emissivity  figures  and  this  fact  could  also  be  used  for  validating  the 
algorithms.  In  this  process,  however,  the  spectrometer  calibration  process  will  have  to  be  evaluated  since  there 
could  still  be  error  of  approximately  0.5K  with  the  calibration  technique  we  are  using.  We  also  conclude  that 
the  surfaces,  especially  the  aluminium  panel  could  be  assumed  to  be  lambertian,  or  at  least  has  reflection 
characteristics  similar  to  these  of  the  infragold  plate  used  to  measure  the  downwelling  irradiance. 

6.0  CONCLUSION 

An  algorithm  for  temperature  and  emissivity  separation  is  proposed  in  this  paper.  It  is  based  on  the  iteration 
on  temperature  principle  where  the  temperature  selection  criterion  is  the  minimum  total  square  error  between 
the  ground-leaving  radiance  and  a  filtered  emissivity  counterpart  of  the  same  radiance.  The  computation  at  a 
given  temperature  involves  the  estimation  of  a  trial  emissivity  that  is  filtered  and  re-used  in  the  computation  to 
determine  a  different  radiance.  The  total  square  error  is  the  sum  of  the  square  of  the  difference  between  the 
measured  radiance  for  each  sensor  bands.  The  main  components  of  the  technique  are  the  filters  and  the  use  of 
the  total  square  error.  The  operational  basis  is  that  when  the  trial  temperature  will  be  very  near  the  sample 
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temperature  the  estimated  emissivity  will  be  smooth  enough  that  the  filtering  will  not  modify  it  appreciably 
and  therefore  the  total  square  error  will  be  minimal  at  that  temperature. 

The  system  is  impeded  by  sensor  noise.  We  obtained  mathematically  an  expression  for  the  estimation  of  the 
impacts  of  noise  in  terms  of  the  standard  deviation  of  the  temperature  estimate.  The  expression  showed  a  very 
good  correspondence  with  Monte-Carlo  simulation  using  the  technique  in  the  same  conditions  as  for  the 
computations.  Experimental  validation  of  the  technique  involves  the  characterisation  of  a  sensor  with  respect 
to  noise  and  the  acquisition  of  a  large  number  of  independent  measurements  and  the  comparison  with  the 
standard  deviation  estimation. 

The  mathematical  model  representing  the  technique's  operations  is  important,  because  it  enables  the 
evaluation  of  the  effect  of  the  main  component  of  the  method,  which  is  the  filter,  with  respect  to  different 
conditions.  The  weather,  the  ground  temperatures  and  the  samples'  emissivities  are  examples  of  these 
conditions  that  affect  the  quality  of  the  results.  It  is  also  important  because  it  could  lead  to  an  approach  for 
filter  design  adapted  to  maximize  the  precision  given  a  particular  condition  of  operation.  Finally,  it  enables  the 
assessment  of  the  result's  quality  once  a  computation  has  been  made  on  a  given  image.  This  particular  aspect 
is  of  interest  for  people  working  in  other  field  of  hyperspectral  image  processing,  because  it  can  conditions  the 
operations  of  detection  algorithms  and  the  subsequent  use  of  imagery. 

A  good  TES  technique  is  also  required  to  lay  the  ground  for  an  autonomous  atmospheric  compensation 
technique,  which  is  the  future  for  high  altitude  thermal  hyperspectral  imagery.  In  this  scenario,  one  possible 
scheme  is  the  use  of  iteration  for  the  atmospheric  profiles  leading  to  credible  emissivity  and  temperature 
estimation  related  to  a  given  pixel's  spectrum. 
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Introduction 


•  What  is  temperature  emissivity  separation 

-  Separation  of  the  ground-emitted  radiance  in  its  two  fundamental  components;  the 
emissivity  and  the  temperature; 

-  Mathematically  ill-defined  problem; 

•  The  unknown  set  always  contains  one  more  element  than  the 
measurement  set; 

-  Use  of  atmospheric  downwelling  irradiance  as  a  reference  for  the  separation. 


Emissivity 

-  Provides  information  related  to  the  nature  of  the  surface 

The  temperature 

-  Provides  information  related  to  the  use  of  the  targets  and  to  their  environmental 
interaction 
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Ru  Radiance  model  (from  8  to  12  microns) 


Transmittance 
Accumulated  path  radiance 


Downwelling  irradiance 
acts  as  a  reference 


Sun  heats  the  surface 
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Processing  chain 


•  Step  1 : 

-  Obtain  from  raw  sensor  data  the  sensor  measured  radiance 

•  Step  2: 

-  Removal  of  path  radiance  and  transmittance  to  obtain  ground 
emitted  radiance 


Step  3: 

-  Separation  of  temperature  and  emissivity 

Step  4: 

-  Process  imagery  to  extract  information 

•  Can  be  done  with  either  radiance  imagery  or  emissivity  and 
temperature  imagery 
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Operation  basics 
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f  Algorithm  structure 


Start  at  a  given 
temperature 


Computes  emissivity 
at  temperature 


Smoothes  the 
emissivity 


Compute  at- 


Pixel  initialization 


Linear  filter 


Temperature 
selection  and 
decision 


Major  concern 

Conditions  the 
algorithm 
operation  and 
performance 
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Radiance  [Watt/(m  .sr.cm  )] 


>  Algorithm  Components 

^  Pixel  Initialization 


Use  2  close  measurements 

-  Assume  same  emissivity 

-  Assume  same  blackbody  radiance  values 


7,  _  L2R1-LIR2 

^  ±  start  - 


(R1-R2Y  (LrL,) 


Constraint: 

Start  temperature  must  be  close  to  true 
temperature. 

High  contrast  between  the  two  chosen 
bands  is  required 
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Algorithm  Components 

Filter 


Most  important  component  of  the  algorithm 

Determines  the  performance 


Performance  depends  on: 

-  Sensor  noise 

-  Atmospheric  radiance 

-  Emissivity  characteristics 

Design  criterion 

-  No  modification  to  original  emissivity 

-  Highly  attenuate  residual  atmospheric  features 
Choose  a  general  low  pass  filtering  system 

Gaussian  shape 
Better  shapes  can  be  found 


Several  lines  of  a  gaussian  filter 
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V  Algorithm  Components 

Temperature  selection 


Minimum  is  found  using 
iteration  on  temperature 
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Temperature  standard  deviation  [K] 


Performance 

Resistance  to  noise 


Performance  depends  on  sensor  noise 

The  true  signal  is  the  contrast  between  the 
ground  radiance  and  the  sky 


The  filter’s  effects  on  noise  influence  the 

temperature  estimation  accuracy.  Narrow 
band  filters  will  be  more  accurate. 


The  lower  the  emissivity  the  lower  the  accuracy 
=>  Temperature  estimation  depends  on  target  self¬ 


emission 


For  a  constant  emissivity;  there  is  no  bias  =>  The  algorithm  behaves  linearly 
introduced  by  noise  for  SNR  above  100 
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Performance 

Other  effects 


•  Variable  emissivity  may  introduce  bias  in  temperature 
estimation 

-  The  filter  must  be  designed  properly 

•  Errors  on  downwelling  irradiance  estimation  introduce  bias 
on  temperature 


-  Radiance  offset  in  incident  light  and  variations 
disruptions  are  responsible  for  that 

-  Filters  properly  designed  may  alleviate  the  problem 

-  Atmospheric  parameters  estimation  must  be  accurate 

Errors  on  atmospheric  compensation  may  well  introduce 
difficulties 
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Measurement  System 


2) 


Methodology: 

Calibration  using  two 
blackbody  radiance  sources 

Sky  irradiance  measurement 
using  infragold  plate 


ABB  Bomem  MR-300 
Aiming  in  a  folding  mirror 
Objective: 

1)  Acquire  ground  truth  data 

2)  Validate  the  algorithm’s 
performance  attributes  in  real 
conditions 
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Extracted  emissity 


Measurement  example 


riJ  Conclusion 


•  We  produced  an  algorithm  for  temperature-emissivity 
separation  for  operation  in  the  Long  wave  infrared  from  8 
to  12  microns 

•  Its  operation  is  based  on  iteration  on  temperature 

•  We  produced  analytical  expressions  for  the  prediction  of 
its  performance  considering: 


The  sensor  noise 
The  environment 

The  emissivity  of  the  targets  under  study 
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Future  work 


•  Establish  experimentally  the  validity  of 
the  analytical  expressions 

•  Validate,  using  real  sensor  data,  the 
algorithm  and  its  performance  prediction 
tools. 

•  Implementation  of  the  algorithm  in  an 
easy  to  use  tool  to  process  data  from  any 
sensor  and  especially  AIRIS 
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